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Abstract 



(N 

> , 

Tlj- \ It is well known that usual quantum teleportation protocols can- 

■ not transport energy. Recently, new protocols called quantum energy 

. teleportation (QET) have been proposed, which transport energy by 

local operations and classical communication with the ground states 
<^ ' of many-body quantum systems. In this paper, we compare two dif- 

. ferent QET protocols for transporting energy with electromagnetic 

field. In the first protocol, a 1/2 spin (a qubit) is coupled with the 
quantum fluctuation in the vacuum state and measured in order to 
^ ■ obtain one-bit information about the fluctuation for the teleporta- 

tion. In the second protocol, a harmonic oscillator is coupled with 
the fluctuation and measured in order to obtain continuous-variable 
information about the fluctuation. In the spin protocol, the amount 
of teleported energy is suppressed by an exponential damping factor 
when the amount of input energy increases. This suppression factor 
becomes power damping in the case of the harmonic oscillator proto- 
col. Therefore, it is concluded that obtaining more information about 
the quantum fluctuation leads to teleporting more energy. This re- 
sult suggests a profound relationship between energy and quantum 
information. 



1 Introduction 



In quantum field theory, the concept of negative energy physics has at- 
tracted considerable attention for a long time. Quantum interference can pro- 
duce various states containing regions of negative energy, although the total 
energy remains nonnegative yj. The concept of negative energy plays im- 
portant roles in many fundamental problems of physics, including traversable 
wormhole |2], cosmic censorship |3], and the second law of thermodynamics 
[1] . In addition, its physical application to quantum optics has been discussed 
[3]. Recently, negative energy physics has yielded a quantum protocol called 
quantum energy teleportation (QET) in which energy can be transported 
using only local operations and classical communication (LOCC) without 
breaking causality and local energy conservation [6]- [8]. QET can be theo- 
retically considered in various many-body quantum systems including l-fl 
dimensional massless fields |H], spin chains [7] and cold trapped ions [S]. 
Based on developing measurement technology with sensitive energy resolu- 
tion for the systems, the QET effect might be observable in future. It may 
be also possible to enhance the amount of teleported energy by preparing 
a large number of parallel QET channels, performing a QET protocol for 
each channel and accumulating each teleported energy so as to achieve de- 
sired amount of total energy. After future experimental verification of QET, 
amazing possibility would be open in principle for nano-technology applica- 
tion of QET. For example, it may be imagined that, without heat generation 
in the intermediate subsystems of the QET channels, energy is transported 
in nano-machines at a speed much faster than the evolution speed of excita- 
tions of the channels. This technology, if possible, helps future development 
of quantum computers in which energy distribution and quantum tasks in 
the devices are completed before heat generation in the system. QET is also 
expected to provide insights on unsolved problems in gravitational physics. 
In fact, a QET process has already been analyzed in black hole physics, and 
from the measured information of zero-point oscillation of quantum fields, 
it can be regarded as controlled black hole evaporation if we consider the 
protocol near the horizon of a large- mass black hole [9] . 

Energy transportation usually requires physical carriers of energy such 
as electric currents and radiation waves. Energy is infused into the gateway 
point of a transport channel connected to a distant exit point. Then, energy 
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carriers of the channel excite and propagate to the exit point. At the exit 
point, energy is extracted from the carriers and harnessed for many purposes. 
On the other hand, in the QET protocols, energy can be extracted from the 
exit point even if no excited energy carriers arrive at the exit point of the 
channel. We locally measure quantum fluctuation around the gateway point 
in the ground state of the channel system and announce the measurement 
result to the distant exit point with zero energy density, where we can ex- 
tract energy from the channel. A key feature is that this measurement result 
includes information about the quantum fluctuation of the channel around 
this distant point via quantum correlation of the ground state of the chan- 
nel system. Therefore, we can infer details about the behavior of a distant 
fluctuation from the result of the local measurement. To compensate the ex- 
traction of this information, some amount of energy must be infused into the 
channel system at the measurement point; this is regarded as input energy to 
the gateway point of the channel. By choosing and performing a proper local 
operation based on the announced information at the distant point, the local 
zero-point oscillation around the distant point can be suppressed relative to 
the ground-state one, yielding a negative energy density. During the opera- 
tion, respecting local energy conservation, positive amount of surplus energy 
is moved from the channel system to external systems. This is regarded as 
output teleported energy from the exit point of the channel. 

One of the important unresolved problems in QET is the theoretical clar- 
ification of the properties in 1+3 dimensions. Protocols in 1-1-1 dimensions 
have already been extensively analyzed in previous studies [6]- [8]. However, 
1+3 dimensional models have not yet been analyzed. In addition, all the pro- 
tocols proposed thus far adopt quantum measurements for discrete-variable 
information. Therefore, it would be interesting to investigate not only a 
protocol with discrete-variable information but also one with continuous- 
variable information. In this study, we carry out a detailed analysis of two 
QET protocols for 1+3 dimensional electromagnetic field in the Coulomb 
gauge. Local measurements of quantum fluctuations in the vacuum state 
of the field require energy infusion to the field. The infused energy is dif- 
fused to spatial infinity at the velocity of light and the state of the field soon 
becomes a local vacuum with zero energy around the measurement area. Ob- 
viously, this escaped energy cannot be taken back to the measurement area 
by local operations around this area if we do not know the measurement 
result of the fluctuation. However, if the measurement result is available, 
we can effectively take back a part of this energy to the measurement area 
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by applying the QET mechanism. By carrying out a local unitary operation 
dependent on the measurement result for the measurement area with zero 
energy density, the fluctuation of zero-point oscillation is squeezed and a neg- 
ative energy density appears around the area, accompanied by the extraction 
of positive energy from the fluctuation to external systems. Needless to say, 
without the measurement result, it is impossible to extract energy from the 
zero-energy fluctuation. One of the two QET protocols we will consider is 
a teleportation in which discrete-variable information about a fluctuation is 
obtained using a measurement with a 1/2 spin (a qubit), and the other is 
a teleportation in which continuous-variable information is obtained using 
a measurement with a harmonic oscillator. The discrete-variable protocol 
is a straightforward extension of the protocol for a 1+1 dimensional field 
proposed in [6]. The measurements are generalized (POVM) ones that use 
probe systems (1/2 spin and harmonic oscillator) strongly interacting with 
local electric field fluctuations during a short time. We prove that for a large 
energy input, the continuous- variable teleportation is more attractive than 
the discrete-variable teleportation. In the discrete-variable case, the amount 
of teleported energy is suppressed by an exponential damping factor when 
the amount of energy infused by the measurement increases. Meanwhile, this 
suppression factor becomes power damping in the continuous- variable case. 
Therefore, it is concluded that obtaining more information about the quan- 
tum fluctuation leads to teleporting more energy. This result suggests a new 
profound relation between energy and quantum information. So far, relation- 
ship between energy and information has been extensively discussed only in 
the context of computation energy cost [10], [11], [12]. The QET viewpoint 
may shed light on a new relationship between amount of teleported energy 
and amount of quantum information about ground-state fluctuations which 
would be characterized by various informational indices including mutual in- 
formation and entanglement. The explicit analysis about this relationship is 
beyond the scope of this paper and will be reported elsewhere. 

The remainder of this paper is organized as follows. In section 2, a brief 
review of the quantization of the electromagnetic field in the Coulomb gauge 
is presented in order to clarify our notations, and the emergence of negative 
energy density is explained. In section 3, we discuss a discrete-variable pro- 
tocol. In section 4, a continuous- variable protocol is analyzed. In section 5, a 
summary and discussions are presented. In this paper, we adopt the natural 
unit c = h = 1. 
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2 Quantization in Coulomb Gauge 



We present a short review of quantization of the electromagnetic field in 
the Coulomb gauge in order to clarify the notations used for later discussions. 
The gauge is defined by 

Ao = 0, divA = 

for the gauge field =(^0; A). Then, the equation of motion of the Heisen- 
berg operator of the gauge field A{t, x) is reduced to the massless Klein- 
Gordon equation given by 

[d^ - V^] A = 0. 

The solution can be expanded in terms of plain-wave modes as follows. 
y(27r)^2|k| fe=i,2 ^ 

where (aj^) is a creation (annihilation) operator of the photon with 

momentum k and polarization h satisfying 

and (k) is a polarization vector satisfying 

Gh (k)* • eh' (k) = 5hh', 

In this study, because we take a sum of two polarization contributions to 
obtain the final results, the reality condition can be imposed on eh (k) for 
simplicity such that ejl^ (k) = e/i(k). In addition, e^j (k) satisfies the com- 
pleteness relation as 

h 

The energy density operator of the field is defined by 
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6(x) = i: (^E(xf +(vxA(x))') 



where E (x) is the electric field operator and :: denotes the standard normal 
order with respect to and . The Hamiltonian is given by the spatial 
integration of the energy density as follows. 



H 



It is a well-known fact that the Hamiltonian (total energy of the field) is a 
nonnegative operator. The vacuum state |0) is the eigenstate with the lowest 
eigenvalue zero of ^ as ^|0) = 0. The expectation value of energy density 
vanishes for the vacuum state as 



|e(x)|0) = 0. 



(1) 



In the later discussion, coherent states are used often. Therefore, we 
present a summary of the related properties of the coherent states. A dis- 
placement operator generating a coherent state from |0) is given by 



U (p, q) = exp 



p(x).A(x)-q(x).E(x) 



(fx 



where p(x) and q(x) are real vector functions satisfying the conditions of 
this gauge as 

V ■ p(x) = 0, 

V • q(x) = 0. 

By using the commutation relation between the gauge field and the electric 
field at time t = given by 

d d ' 

i\5ab- ) 5 (x - y) , 



(2) 



i„(x), h{y) 

it is easily verified that IJ (p, q) displaces A(x) and E(x) as 
m (p, q) E(x)t/ (p, q) = E(x) + p(x), 



m (p, q) A(x);7 (p, q) = A(x) + q(x). 



(3) 
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In addition, we are able to prove a product formula as 
U (Pi,qi) U (P2,q2) 



exp 



(Prq2-qi-P2)c?^3; 



U (pi + P2,qi + q2) • 



This implies that the set of U (p, q) forms a unitary ray representation of 
the displacement group of the field. Coherent states generated by U (p, q) 
are defined by 

|(p,q)) = f/(p,q)|0). (4) 
By using the Fourier transformation of p(x) and q(x) defined by 



P (k) = / p(x)e-^''-^d3x, 



Q (k) = j q(x)e-''^-"d3x, 



the coherent states are explicitly written in terms of the creation operator 
as follows. 



(p,q)) 



= exp 



X exp 



(27r)'2|k| 
^J{2^^f2\\^\ h 



|P(k) -i|k|Q(k) 



J]e,(k)-(P(k)-i|k|Q(k))a; 



|0). 



Prom the above expression, it is easy to prove that | (p, q)) is an eigenstate 
of the annihilation operator d^ such that 



«kl(p,q)) = 



=e,(k).(P(k)-i|k|Q(k))|(p,q)). (5) 



V(2^)'2|k| 

The inner product of two coherent states is exphcitly calculated as 
((Pi,qi) I (P2,q2)) 



: /(pi-q2-qi-P2)<i^^ 



X (i^^p|^|Pi(k)-P2(k)-i|k|(Qi(k)-Q2(k))|2 



(6) 
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Next, we examine the emergence of a region with negative energy density 
in this standard theory. As a simple example let us consider a super- 
position state l^') of the vacuum state |0) and a two-photon state |2) such 
that 



1^) = cos^lO) +e^''sin^|2), 

where 9 and 5 are real parameters with 6 G [0, vr] and 5 G [0, 27r]. Generally, 
an off-diagonal element of the energy density (0|e(x)|2) does not vanish 
for a fixed point x. This is because e (x) includes a non- vanishing term 
proportional to a^al^i- This fact indicates the emergence of negativity of the 
energy density as follows. The expectation value of the energy density for 
the state |\E') is calculated as 



(^|£(x)|v&) 

= 2cos^sin6'(cos(5Re(0|£(x) |2) - sin(5Im(0|£ (x) |2)) 
-fsin2^(2|e(x) |2). 

In this result, let us set the parameters 6 and 5 so as to satisfy 

(2|e(x)|2) 



cos 9 



sin 9 



tan 5 



(2|£(x)|2)2 + 4|(0|5(x)|2)|^ 
2|(0|e(x)|2)| 

(2|£(x)|2)2 + 4|(0|e(x)|2)|^^ 
lm(0|£(x) |2) 



Re(0|£(x) |2)' 

Then, {^\e (x) |\E') is evaluated as a negative value as follows. 



(^|e(x) 1^) 



(2|5(x)|2)2 + 4|(0|£(x)|2)|^-(2|e(x)|2) 



< 0. 



Therefore, the emergence of negative-energy regions is not a peculiar phe- 
nomenon in quantum field theory. Quantum interference in the superposi- 
tion of photon number eigenstates yields negative values. It should be re- 
emphasized that despite the existence of regions with negative energy density. 
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the expectation values of H remain nonnegative. This imphes that there ex- 
ist regions with a sufficient amount of positive energy so as to make the total 
energy nonnegative. As described in sections 3 and 4, this negative energy 
plays a crucial role in the QET protocols. 



3 Discrete- Variable Teleportation 

Our protocol for QET with a 1/2 spin probe is a straightforward extension 
of that in ^ and comprises the following three steps. 

(1) At time t = 0, the spin probe is strongly coupled with the vacuum 
fluctuation of the electric field within a finite compact region Vm during a 
very short time. In this process, some information about the fluctuation is 
imprinted into the probe. Positive energy is inevitably infused into the field 
during the measurement process, as seen later. The amount of energy is 
denoted by Em. When this energy is infused, positive-energy wave packets 
of the field are generated and these propagate to spatial infinity with the 
velocity of light. 

(2) After switching off the interaction, projective measurement of the z- 
component of the spin is carried out. If the up or down state is observed, 
we assign s = + or — , respectively, to the measurement result. This implies 
that we obtains one-bit information about the field fluctuation via the probe 
measurement. 

(3) At time t = T, it is assumed that the measurement has finished and 
the wave packets have already escaped from the region. Hence, the energy 
density in the region Vm is exactly zero. Then, a local displacement operation 
is carried out depending on s within Vm- Even though we have zero energy 
in Vm, positive energy is extracted from the field fluctuation during the local 
operation, generating negative-energy wave packets of the field. The amount 
of negative energy of the wave packets is denoted by Eo{= — \ Eo\). Therefore, 
the amount of energy extracted energy from the field is given by + lE'ol- 

In step (1), the measurement interaction between the electric field and 
the spin probe is given by 
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H^it) = g{t)a,G, (7) 

where is the 2;-component of the Pauh matrices of the spin probe; g{t), a 
time-dependent real couphng constant; and G, a Hermitian operator defined 
by 
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G = - - / a„(x) ■ F.{^)d?x. 



(8) 



Here, am(x) is a real vector function with a support equal to Vm satisfying 
V • am(x) = 0. In addition, by taking a short-time limit for switching the 
interaction, we set 



g{t) = 5{t-Q), 



(9) 



The initial state of the spin probe is assumed to be the up state of the 
x-component of the spin given by 



In step (2), the measurement operators M± [T3] are defined by 



M± = (±,|Texp 



'o 



(10) 



where \—x) is the down state of the x-component of the spin given by 



—I 

71 



1 

-1 



By using Eq. ([9]), M± are computed as 



M+ = cos G, 
M_ = sin G. 



(11) 
(12) 



Using Eq. (j4]), the post-measurement states of the field obtaining s = ± are 
calculated as 
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|^+) = -=M+|0) = —— [e^^ (0,aj) + e~^^ (0, -a, 



1^-) = ^Af_|0) = [e^f I (0,aj ) - e-f I (0, - a„))] , 

where j9± is the measurement probabihty of |±^) for the spin probe and it is 
evaluated as 

p± = (0|MiM±|0) = \. 
After the measurement, the average state of the field is given by 



The amount of energy of the field after the measurement is not zero but 
is instead a positive value given by 

Em = y Tr [pMsi^)] rf^x = ^ y (V X a„(x))' £x. (13) 

This evaluation of Eq. f|T3l) is achieved by using Eq. ([5]). Because the initial 
state of the field is the vacuum state and it has zero energy, the positive 
value in Eq. f lT^ implies that the manipulation of the measurement requires 
positive work from outside. In general, any local operation on the vacuum 
state infuses nonzero energy into the field if the post-operation state is not the 
vacuum state, because the Hamiltonian is a nonnegative operator. In Figure 
1, the measurement process of QET (step (1)) is depicted in the xy-plane 
slice. The measurement area Vm is represented by a circle in the plane. 

Next, let us discuss time evolution of the field after the measurement. 
The average state evolves as 



Pm(t) = J]p.^7(t)|^«)(^.|^7(t)^ 



by the time-evolution operator U{t) = exp 



-itH 



. The Heisenberg opera- 



tors A(t, x) and E(t, x) are given by the Schrodinger operators A(x) ( = A(0, x 
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and E(x) ^= E(0,x)j as follows. 

A(t,x) = J A(2)(t,x-y)A(0,y)d3?/ 

+ j A(i)(t,x-y)E(0,y)A, (M) 

E(t,x) = J dA^'\t,^-y)A{0,y)d'y 

+ J A(2)(t,x-y)E(0,y)d3y, (15) 

where A*^^^ and A*^^) are Lorentz-invariant functions defined by 



A^'nt,x-y) = i-e(t)5(t2-|x-y|^), 

A(^)(t,x-y) = a,A«(t,x-y). (16) 

Both the functions A*^^^ and A*^^) vanish in the off-light-cone region with 
t^ — |x — yl^T^O and satisfy the massless Klein-Gordon equation: 

[d^-V'] A(")(t,x-y) = 0. 

Substituting Eq. (HM and Eq. (|T5l) into i{t, x) yields the time evolution of 
the average value of the energy density as 

{i{t, x)) = Tr [pUm^)] = \ [n(t, x)' + b(t, x)'] , (17) 
where n(t,x) and b(t,x) are given by 

b(t,x)= j A(2)(t,x-y)Vxa„(0,y)t/32/, (18) 
n(t,x)= j 9iA(2)(t,x-y)a„(0,y)A- (19) 

Taking account of the explicit form of A*^^^(t,x — y) in Eq. f|T6|) . the above 
result teaches us that the wave packets excited by the measurement soon 
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escape from the measurement area to spatial infinity with the velocity of 
light. This ensures that the energy density in Vm returns to zero. At time 
t = T, the state PmiT) is assumed to be a local vacuum state with zero 
energy density: 

Tr[p„(T)e(x)] = 0, 

for xGV^- 

In step (3), let us consider a local displacement operation within Vm 
depending on s defined by 



tJs = exp 



is9 / fo(x) • A(x)(i^x 



(20) 



where ^ is a real constant fixed below, and fo(x) is a real vector function with 
the support Vm that satisfies 

V ■ fo(x) = 

so as to preserve gauge invariance of (Jg. In this process, positive energy is 
released on average from the field to the apparatus executing Ug by taking a 
proper value of 6. Let us introduce the localized energy operator around Vm 
as 

Ho = w{x.)i (x) d^x, 



where w(x) is a real window function of Vm that satisfies w(x) = 1 for xGKn 
and decays rapidly outside Vm- The average state after the displacement 
operation is given by 

p = ^ UgUiT)Mg\0){0\MlU^iT)Ul 

s=± 

For this state, the average energy of the field around Vm is defined by 



^n = Tr 



Hop] = Y,{0\MlU\T)UlHoUsU{T)MM- (21) 



Here, the operator tj\T)UlHotJs(j{T) is rewritten using Eq. ^ and Eq. 

as 
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U^{T)UlHoUsU{T) 
w{x)i (T, x) cPx 



1 



+ s9 J f,(x) ■ E(r,x)rf3x + -9^ J i,{x.y£x, 

where e(r, x) is the Heisenberg operator given by U'^ (T)i (x.) U(T), and we 
have used w(x)fo(x) = f^lx). Substituting the above relation into Eq. fET]) 
yields the following relation. 



s=± I--' 



X 



=± 



fo(x) ■E(T,x)c/3a; 
J^MlMs] |0). 



Ms\0) 
M,|0) 



(22) 



In order to further simplify the form of Eq, the following equation is used. 



f„(x) ■ E{T,x)d^x, j a„(y) • E{0,y)d'y 
J J Ux)-Si,^iy)dTA^^\T,x-y)d'xd\ 



(23) 



This is proven by Eq. (fT5|) . Taking account of the supports of fo(x) and 
a„i(y), it is verified that the integration of the right-hand side in Eq. (123|) 
is zero because Ot^^'^^T, x — y) vanishes when — |x — y|^ ^ 0. Thus, we 
obtain the relation 



fo(x) ■ E(r, x)d^x, / a„(y) ■ E(0, y)d^y 



0. 



(24) 



In a similar manner, we can show the following relation using Eq. (IT^ and 
Eq. m- 



w(x)£ (T, x) d^x, / a„(y) ■ E(0, y)d^y 



0. 



(25) 
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From these relations, we can show the commutation relations given by 



w(x)£(T,x)rf^x, Ml 



fo(x) ■E(r,x)(/3x, Ml 







by recalling that M] is a function of / am(y) ■ E(0, y)d^y, as seen in Eq. (|8]), 
Eq. ( II ip . and Eq. (I12p . Therefore, we can move the positions of Mj in the 
first and second terms of the right-hand side of Eq. ( 12^ to the right, and 
rewrite Eo as 



Eo = (0| j w(x)£ (T, x) |0) 
+ ^(0| (^j fo(x)-E(T,x)rf3z^ (^sMlMs^ |0) 

By using two relations of the measurement operators given by 
J2 Ml Ms = cos^ G + sin^ G = 1, 

s=± 

J2 sMlMs = cos^ G - sin^ G = cos (^2G^ , 

s=± 

the following expression of E^ is obtained. 

Eo = (0| j w(x)£ (T, x) c?3x|0) + er]+ ^9'^, 



where rj and ^ are real constants defined by 



V= {0\{ / fo(x) ■E(T,x)(i3x cos(2G) |0), 



(26) 



(27) 

(28) 
(29) 
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The first term of tlie riglit-liand side of Eq. fl27l) vanishes because we can 
transform it into 



j w(x)£(T, x) d^x\0) 
{0\U\T) j w(x)£(0,x) £xU{T)\0) 
w(x)(0|£(x) \0)d^x = 0. 



Here, we have used f/(T)|0) = |0) and Eq. ([T]). The expression of 77 in Eq. 
can be further simplified by noting that 



(30) 



cos (2G) |0) = - [I (0,2a^)) - | (0, - 2a^))] . 
By using Eq. ([5]) and Eq. ( l30l) . the following relation is derived. 

(o|E(r,x)cos (20^ |0) 

= -(0|(0,2a^)) j d^THT,^-y)siJy)£y, 
where A {t, x) is a Lorentz invariant function defined by 



(31) 



A(t,x) = J 



(27r)^|k| 



cos (k • X — |k|t) . 



Contrary to A*^^) and A*^^^ (t.x), the function A (t,x) does not vanish 

in the off-light-cone region: 



A(t,x) 



1 1 



|t2-xV0 



27r2 t2 _ x2 

The factor (0| (0,2am)) is real and calculated from Eq. (E]) as 



(0|(0,2aJ)=exp 



where a^ (k) = / am(x)e ""^'^d^x. Substituting Eq. (|3T]) into Eq. fl28|) gives 
the final expression of 77 as 




0,2a™)) / / a^A(T, X - y)f„(x) ■ a.,r.{y)d^xd^y. (32) 
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Minimization of Eo with respect to 6 is achieved by taking the parameter 9 
as 

9 = -^. (33) 

Substituting the value of 9 into Eq. ( 1271) yields the following expression of 
Eo. 

E. - (34) 

Substituting Eq. f l2I?]) and Eq. fl5^ into Eq. flM|) . we obtain our main result 
in this section: 

^ ^ _^ [//a|A(r,x-y)f»-a^(y)rf3xd3^]^ 
where Dq is an exponential damping factor with respect to such that 



|(0|(0,2a^))|' = exp 



|^|a„(k) 



(36) 



Because ^ > 0, as seen in Eq. ( 129|) . it is of importance that Eo is negative in 
Eq. (El: 



= - l^ol < 0. 

Respecting local energy conservation, this result implies that positive en- 
ergy + 1 -Eo I moves from the field to external systems during the displacement 
operation Us because the energy of the field around Vm is zero before the 
operation. Nonnegativity of the total energy of the field ensures that \Eo\ is 
smaller than E^, as discussed in reference [6]. The wave packets with nega- 
tive energy begin to chase after the positive-energy wave packets generated 
by the first measurement. Figure 2 shows the energy extraction process (step 
(2)). The outside doughnut- shaped region represents the propagating wave 
packet generated by the measurement in step (1). The displacement opera- 
tion Us is performed in the measurement area Vm denoted by the circle in the 
figure. In Figure 3, the average energy density is plotted immediately before 
the displacement operation at t = T — 0. We have zero energy around the 
area where Us is performed. In Figure 4, the average energy density after the 
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operation is plotted. Negative energy density appears around V^. In conclu- 
sion, a part of the escaped energy is effectively retrieved to the measurement 
area by this QET protocol. 

In this discrete- variable case with the spin probe, the amount of extracted 
energy \Eo\ in Eq. fl35l) is suppressed by the exponential damping factor Dg in 
Eq. fl5B]) when the magnitude of and the infused energy Em increase. This 
suppression becomes power damping for a continuous-variable teleportation 
with a probe harmonic oscillator, as seen in the next section. 



4 Continuous- Variable Teleportation 



In this section, we analyze a QET protocol with a continuous variable. 
The essential part of the protocol is the same as that discussed in section 3 
except the probe system and the operation dependent on the measurement 
result. As the probe to measure the quantum fluctuation of electric field, let 
us consider a harmonic oscillator. The free Hamiltonian is given by 



where q is the position operator and p, the momentum operator. This os- 
cillator couples with the electric field by the measurement interaction such 
that 



operator in Eq. (|8]). The initial state of the oscillator is the ground state 
\g) of Hho- Introducing an eigenstate \q) of q with eigenvalue q, the ground 
state is described by 



After the measurement interaction is switched off, we measure the position 
q. The measurement operator for measurement result q is computed as 



Hho = ^P' + 2f 



(37) 





1/4 



exp . 



Mg = {q\ exp -ipG \g). 
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More explicitly, Af„ is written as 



(38) 



The operator Mq satisfies the following relations. 

/oo 
MiMqdq = 1, 
'OO 

/oo 
qMlMqdq = G, 
-oo 

/°° ^ ^ - 1 

q'MlMqdq = G' + -. 
-oo ^ 

The state after the measurement is given by 

Mq\0){0\Mldq. 



(39) 
(40) 
(41) 



In this measurement, energy is infused into the field as given by the protocol 
described in section 3. The expectation value of the infused energy density 
is evaluated as 



(£(t,x)) = Tr [p:„£(t,x)] = / {0\Mli{t,^)Mq\0)dq 



(42) 



where e(t,x) = W{t)e {:>i)U{t). The value (£(t,x)) can be calculated by 
using the following integral formula. 



exp (-^2) 



1 



e 4 



'tpq 



dp. 



(43) 



From Eq. (H3!l . the following relation is derived. 

1/4 



M,|0) 



exp 



2V/' 1 



-{i-gY 



dpe~ 



|0) 



exp 



ipG 



|0). 



From Eq. Mq\0) is rewritten as 
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2\^/^ 1 



M,|0) = ( - 1 -= / dpe---*n«-4j| (0,pa^)). (44) 



Substituting Eq. (jUj) into Eq. yields the following manipulation: 



(e(t,x)) 
1 



V27r 
1 

1 

((0,a. 




e-^-^ ((0,pa„) |£(t,x)| (0,pa^)) (^j e'(P-P)ii-^) dpdp 
(ipe 2 ((0,pa^) |£(t,x)| (0,pa^)) 



p^e 2 dp((0,a„) |e(t,x)| (0,a^)) 
|e(t,x)| (0,a„)), 



(45) 



where we have used {{0,pam ) |£(t,x)| (0,pa^)) =p2((o,a„) |£(t,x)| (0,a„)). 
From Eq. fH5l) . we can show that the input energy (£(0,x)) d^x in the 
measurement process is the same as that of the discrete-variable case de- 
scribed by Eq. f[T5]) of section 3. Besides, in the same manner as the deriva- 
tion of Eq. (fT7|) . it is possible from Eq. (H5|) to show the equation 

(£(t,x)) = l[n(t,x)Vb(t,x)^], 

where b(t,x) and n(t,x) are given by Eq. f[T5]) and Eq. flTI?]) . respectively. 

After the wave packets generated by the measurement escape from the 
measurement area, we execute at t = T, when the field has zero energy inside 
Vm, a local unitary operation dependent of the measurement result q given 
by 



exp 



iqO' / fo(x) ■ A{x)£x 



(46) 



which is analogous to Eq. ( |20|) . The parameter 9' is real and fixed later. 
Then, the energy around the measurement area is evaluated as 

E'^= f dq{0\MlU\T)UlHM,U{T)M,\0). 



By using Eq. ([39])-Eq. gl]) and 
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U\T)UlHoU,U{T) 
= J w(x)£(T,x) d^x 



1 



e'q j fo(x) • E{T,x)£x + ^9'\^ j U^fd^x, 

the average energy E'^ is rewritten, in a manner similar to that in section 3, 

as 

E'o = e'v' + lo"{o\ {g' + \^ |o), 

where ^ is given by Eq. ( l29i) and t]' is a real constant that is evaluated as 



fo(x) ■ E(T,x)d=^x G|0) 



-(0| i^j f„(x) ■ E(r,x)d3xJ a„(y) • E(y)d3yJ |0) 

i y y a^A(T, X - y)f„(x) ■ i,miy)d'xd'y. (47) 



By fixing the parameter 6' as 



the minimization of E'^ with respect to 9' is attained as 

^'2 



e(0|(G2 + i)|0) 



(48) 



E' 



By substituting Eq. ( l29l) and Eq. ( H71) into the above equation, we obtain 
the final expression of E'^ as 



Eg = —Dho 
where Dho is given by 



[/ / 91 A(T, X - y)f „(x) ■ ^,,iy)d^xd^y] ^ 



2/fo(x)^rf3 



(49) 



D 



ho 



4^2 + 1) |0) 



l + - + 2y|k||a„,(k)r^ 



-1 



• (50) 
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It is noticed that the form of E'^ in Eq. fl49p is the same as that of Eo in 
Eq. fl35p . except for the form of damping factor Dg. Unhke the exponential 
damping of Dg, the suppression factor Dho is a power damping factor with 
respect to |am(k)|. Because of this weak damping, the extracted energy 
\E'J by this protocol does not vanish even for a large amplitude limit with 
IsLm fk) I — 7- oo as 



\E'\ 



/|ki|cos(|ki|r)n^(kO*-fo(kO0 



4 / 



d^k2 



where fo (k) = J fo(x)e ^^'^cPx and (k) is a rescaled amphtude given by 



(k) 



(k) 



/ |k'| (k'^ 



1 2 d^k' 



As a conclusion, we can say that obtaining more information of the fluctu- 
ation leads to teleporting more energy. It is an important question what 
measurements and operations attain the maximum transporation rate of en- 
ergy in the QET mechanism, however, remains unsolved yet. 



5 Summary and Discussion 



We have analyzed in detail two protocols of QET for the electromagnetic 
field and shown that a part of the lost energy in the measurement can be 
retrieved by use of the measurement result. The amount of energy infused by 
the measurement is the same in both the cases and it is given by Eq. (fT3|l . 
For the discrete-variable case, the measurement with the operator Mg given 
by Eq. ( iTTl) and Eq. ( IT2l) is performed. The amount of retrieval energy \Eo\ 
is given by Eq. ( I35l) for the displacement operator in Eq. ( I20l) with 9 fixed 
in Eq. f l33|) . For the continuous- variable case, the measurement with the 
operator given by Eq. f p8|) is performed. The amount of energy retrieved 
\E'J is given by Eq. f l49|) for the displacement operator in Eq. fH6|) with 9 
fixed in Eq. (148|) . For large amplitude la^ (k)| and large energy input Em, 
the continuous-variable teleportation is found to be more preferable than 
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the discrete-variable teleportation. In the discrete-variable case, the amount 
of extracted energy is suppressed by the exponential damping factor Dg in 
Eq. fl36l) when the energy infused by the measurement increases. For the 
continuous-variable case, the suppression factor becomes the power damping 
factor Dho in Eq. f l50|) . Therefore, it is concluded that obtaining more 
information about the quantum fluctuation leads to teleporting more energy. 

In future QET experiments, the separation between Vm and the region 
of the escaping wave packet with positive energy should not be so large to 
extract an observable amount of energy. To see this, let us take a large- 
separation limit as T ^ |x — y| for x, y G V^. Then, the amount of tele- 
ported energy decays as 

l^oloc^, 

in both the protocols. This rapid decay becomes one of the serious obsta- 
cles to observing the extraction of energy by QET. Therefore, for the best 
implementation of the QET in 1+3 dimensions, T should be of the same 
order as the measurement area size. From this viewpoint, QET protocols are 
more attractive for physical systems described effectively by H-l dimensional 
massless field models that have slower-decay properties {Eo oc T"^) of the 
transported energy for a large separation T, as discussed in references [6]- [7]. 
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Figure Caption 

Figure 1: The first schematic diagram of QET in the xy-plane shce. The 
local measurement is performed with the infusion of positive energy to 
the field and the measurement result s is obtained. The measurement area 
Vm is represented by a circle in the plane. A positive-energy wave packet is 
generated in the field system and it escapes to spatial infinity at the velocity 
of light. 

Figure 2: The second schematic diagram of QET in the xy-plane slice. 
The outside doughnut-shaped region represents the propagating wave packet 
generated by the measurement. By using the measurement result s, the dis- 
placement operation Ug is carried out in the measurement area Vm denoted by 
the circle. After the operation, a wave packet with negative energy — |i?o| is 
generated, accompanying the extraction of positive energy + |£^o| to external 
systems. 

Figure 3: The average energy density is schematically plotted immedi- 
ately before the displacement operation at t = T — 0. The positive-energy 
wave packet generated by the measurement has already escaped from the 
measurement area where Us is performed. 

Figure 4: The average energy density is schematically plotted after the 
operation. The wave packet with negative energy — |-Eo| is generated by the 
extraction of positive energy + |-Eo| to external systems. 
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